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Abstract 

In this paper we consider classes of nonlinear systems for idiich the optimal 
(minimum variance) estimator is finite dimensional. Finite dimensional optimal 
nonlinear state estimators are derived for bilinear systems evolving on nilpo- 
tent and solvable Lie groups. These results are extended to other classes of 
systems involving polynomial nonlinearities. Finally, the concepts of exact 
differentials and path-independent integrals axe used to derive optimal finite 
dimensional estimators for a further class of nonlinear systems. 



1. INTRODUCTION 

It is well knovm that the class of linear dynami- 
cal systems with linear observations and white 
Gaussian plant and observation noises is particu- 
larly appealing, because the optimal state estima- 
tor consists of a finite dimensional linear systen 
(which is easily implemented with the aid of a 
digital con^juter) . In general, the optimal (mini- 
mum variance) estimator for a finite dimensional 
nonlinear system consists of an infinite dimensio- 
nal systen of moment equations, and approximations 
must be made for practical implementation. 
Consequently, one is led to inwstigate subclasses 
of nonlinear systems which admit finite dimensional 
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Optimal estLnators. Such an investigation not 
only identifies syst«»)s for which optimal estima- 
tion is computationally feasible, but it also pro- 
vides valuai.'le theoretical insight into the under- 
lying structure of optimal estimation for general 
nonlinear systems. 

There is, in fact, a class of nonlinear systems 
>idiich possesses a great deal of structure — the 
class of bl) inear systems. Such systems have been 
studied by several authors t end many tools 

from the theories of Lie groups and differential 
geometry have proved to be- quite useful. Estima- 
tion for bilinear systems on abelian Lie groups is 
discussed by Lo and Nillsky [6] f and some optimal 
finite dimensional estimators are derived; these 
results are generalized to a larger class of 
systems by Willsky (91 . 

In this paper wo considcir bilinear systems evoJving 
on solvable arwl Lie groups. fat: but;li 

systems, we — H J trrrrv— ' b.t mjeirri''’ 'f.nrJ ’•ti— «» 
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stochastic systssi driven by the innovations. 

These results are extended to other systems with 
polynomial nonlinearities. In addition/ the con- 
cepts of path-indepeiulent integrals and exact 
differential forms are used to prove the exis- 
tence of finite dimensional estimators for ano- 
ther class of nonlinear systems. 

2. HOTATION AND FOmJLATlON OF THE PROBLEM 

The basic bilinear equation considered in this 
paper. is 

H 

ict) ° (A. + I U. (t)A )X(t); X(0) « I (1) 

'■ i-1 ^ ^ 

vdiere the A. are given k x k matrices* X is a 
k X k matrix, is the i con^nenent of u, and 
u is the n-dimensional "colored noise” process 
generated by the finite dimensional linear 
system 

dx(t) • P(t)x(t)dt + (t)dw(t) (2) 

u(t) - C(t)x(t) O) 


Here w is a standard Brownian motion process, 

Q ^ 0, x(0) has a Gaussian distribution indepen- 
dent of w, and the pair [F(t),, is 

stabillzable [IB]. The existence of a unl^e 
•olution to (l),(2) is proved in [15]rU6]> 
Although X by itself is jiot a Markov process, it 
can be shown that the pair (X,x) is Markov. 

As in the deterministic case [i] , the solutior X 
of (1) evolves on a matrix Lie group. More spe- 
cifically, we define to be the mati ix 

Lie algebra generated by {a^, i=0, 1 ,*..,n}; i.e.. 


SB is the smallest subspace of k x k matrices con- 
taining {a^, i*=0, and closed under tha 

commutor product lP,Ql = PQ-QP. We also defina 
the matrix Lie group G ® {exp^}_ associated with 
^to be the smallest group (under matrix' multi- 
plication) containing exp L for all L ei?. If 
X(0) C G, then X (t) EG for all t ^.0. 


In the sequel we will be primarily concerned 
with systems in which SB is & solvable or nilpo - 
tent Lie algebraj such systems evolve on solvable 
or nilpotent Lie groups. 

DefinitioiT ! ♦ A Lie algebra SB is solvable if 
the derived scries of ideals 1231 




OJ? the 
POOB 


n ^ tezBlnates in {o}. SB is nilpotent if the 
lower central series of ideals 


JP” - ^ {[L^,L2lb^e^. > 1 

texsiinates in {o}. SB is abelian if ^ “ {o}. 

We state here two results conseming canonical 
representations of nilpotent and solvable Lie 
algebras vrtiich will be particularly germane to our 
study: the reader is referred to (231 for further 
, propexrties of solvable and nilpotent Lie algebras 
and groups. I.et ^denote the con 5 >lex numbers, 
gt(n, denote the space Of n x n matrices with 
coB^lex entries, and GL(n,^) denote the space of 
nonsingular complex n x a natrices . 

Lwaaa 1 (23 , p.214] Let ^be a Lie algebra 
of matrices in gi(n,^) . Then iPis solvable if 
and only if tl^re exists a matrix P e GL(n,Sf). 
such that. £o: all elements A C SB, 'the ciatelx 
S - PAp"^ is n upper triangular' form 
» 0 for i. > j) . 

Lenwa 2 12 t, p.224]: Let i? be a Lie algebra of 

matrices in gi(n,^). Then JPis nilpotent if and 
only if there exists a matrix R E GL(n ,^) such • 
that, for all el^nts A e^, the matrix B » PAp"^ 
has the block diagonal . form 

^j^(A) * 

0 •. 

4>, (A)J . 

* 



The functions S are linear. Furthermore, 

=0. 

The block form for nilpotent Lie algebras will be 



called the nili>otent canonical forn . 

The first estimation problem to be discussed In** 
volves the state equations (1)~<3) and the p— di- 
mensional observation process 

d 2 {t) » H(t)x(t)dt t (t)dv(t) (4) 

where R > 0 azKi v is a standard Brownian notion 
independent of w and x{0) . This observation pro- 
cess is of interest in the problem of estimating 
the attitude of a rotating rigid body by titeans of 
a strapdown inertial navigation system [5) *173. 

The criterion for the optimal estimate 
(X(t{t) ,x(t|t) } will be the minimization of the 
conditional error covariance 

B[<x(t) - x(t|t>) ’ (x(t) - x(tlt)) 

■f tr{(X(t)-X(t|t)J ’ (X(t) -X(t|t))}| 

where "tr" denotes trace and « {* (s) ,0<s<t}. 

It is well 3cnown [14] that the causal mininura- 
varlance estimates are given by the conditional 
means 

X(t]t) » E*^[X(t)l ~ E[x(t) |z^I 

i(t|t) «■ E^[x(t)3 - EUU)!*^] 

(wa will use the three notations for conditional 
expectation interchangeably) . The computation of 
)c(tjt) is performed by the finite-dimensional 
(linear) Kalman, filter; moreover » the conditional 
density of x(t) given is Gaussian with mean 
x(t|t) and nonrandom covariance P(t) [14]. As 

A . 

remarked in Section I, the computation of X(t|t) 
requires in general an infinite dimensional systeo . 
of equations. We will show in the succeeding 

f 

sections that X{t|t) can be computed with a 
finite-dimensional nonlinear estimator if (l)-(3) 
evolves on certain nilpotent or solvable Lie 
groups. 

3. REDUCTIOn OF THE GEIXERRL PROBl^ 

In this section we show t lat some estimation pro-^ 
blems on solvable Lie groups can be solved by 
considering an estimation problem on a particular 
nilpotent Lie group. The first lemma generalizes 
a result of VJilisky I'j] ■ 'fh'-- *s . 

f',r i «, f.'i' _ ‘-' 
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9A(kr9f)by 

ad,(B) ■ AB-£A - [A,8l <5) 

A 

The notation adj' denotes the i^*' power of the 
A 

operator ad^. 

I^enna 3t consider the equations (l}-(4), and 
letJSf^ be the ideal in spanned by 

ad^ (A.) j-=L,...,N; i»0,...,k^-l. 

*0 3 

Define the k x k matrix valued process 
-A-t 

Y(t) • e ° X(t) (6) 

Then there exists a matrix D(t) such that Y satis- 




V * 


is a basis for and 


y(t) « D(t)x(t) 


In addition, X can be coaqnited according to 


X(t|t) 


e ° Y(tlt) 


Leimaa 3 enables us, without loss of generality, 
to examine the estimation problem for Y (t) 
evolving on the normal subgroup «{ exp.^^}^, 

rather than for X(t) evolving on the full Lie 

group 6 •> {exp^}_. The particular case with which 

we will be concerned is that in which ^ is solvable 

and is nilpotent. In fact, it can easily be 

shown that if is nilpotent, then i? must be 

solvable I how' ver, the converse is not true. 

According to ..emma 3, for such systems we need only 

consider the case in which A « 0 and iP. is 

o o 

nilpotent. 

Exawg)le 1 : Assume that A^ is upper triangular, 

amd {a^,...,A^} are strictly upper triangular 

(diagonal elements are zero) . Then SC is solvable 
and is nilpotent. 


Cxatnplo 2: Tn this Plnitn cllmonn 1 fi.ii ot 



By ntsdns of Loireoa 2, the problem can be further 
reduced to the consideration of Lie algebras In 
nllpotent canonical form. 

I*PPa 4; Consider (l)-(4), where A ■ 0 

. o 

and S is nilpotent. Then there exists a matrix 
P e GL(k,^) such that 

X(t|tJ » PY(t[t> 

where Y satisfies (7) and are in 

nllpotent canonical form. 

l>y means of the following trivial lemma, 
we reduce the problem to the consideration of one 
block in the nilpotent canonical form. 

Lemma 5 ; Consider (l)-(4), where A » 0 airi 

o 

. ,Ajj} are in nllpotent canonical form. 

Then X{t) has a block diagonal form conformable 
with that of {Aj^,...,Aj^}. 

Thus the system (1) can be viewed as the direct 
sum of a number of subsystems; for each 
k^-dimensional subsystem, e gn(kj^) 

(here we have defined gn(m) to be the Lie 
subalgebra of gil(in,^) consisting of the upper 
triangular matrices with egual diagonal elements) . 
He now state the major theorem on finite diiaen- 
slonal estimation for such 'systcwos. 

Theorem 1 ; Consider (1)~(4), where A « 0 

o 

and {a^,...,a^} c gn(k).. Then the conditional 

mean X(t|t) can be confuted by a finite-dinensio- 
nal system of nonlinear stochastic differential 
equations. 

Theorem 1 is proved by induction; the case k=3 
will be proved in Section IV, and the induction 
step is proved in (16J . We note that our result 
also includes the result in [8] as a special case. 

4. PROOF FOR k»3 

For simplicity of notation, the theorem will be 
proved for the case C(t) = I; the proof is pre- 
cisely the same for arbitrary C(t). For k-3, 

we assume that n«N'°4, and A « 0, A, a i, 

o I 


0 

1 0 


0 0 l" 


0 0 o‘ 

0 

0 ( 

A - 

0 0 0 


0 0 1 

0 

0 { 


; “ “J 

0. 0 0 


Then {a^, is a basis for gn(3). The 

solution of (1) can be expressed! in closed form as 


X(t> 


«diere 


y.{t) y (t) I 

» e y,(t) e (y,(t)+C(t))| 


y, (t) 

0 e « (t) 


y^ct) 


z 

- I 


(s)ds 




,t!i 


C(t) 


0 0 


m 


(11) 




Renaricj In the sequel we will apply, without 
further comment, a version of the Pubini theorem 
117] ^Ich permits the interchange of integration 
and conditional expectation, , Since we are'dealing 
only with, integrals of .products of Gaussian random 
processes, th< use of the Pubini theorem is easily 
justified [16 . 

It is evident from (10) that the computation of 
X(t|t) (in pa-'ticular, X^^f^lt)) requires the 
knowledge of i.he statistics of the entire process 
{x(s) , 0 £ s < t} given 2 . Therefore we first 
define the conditional cross cross-covariance 

P(0^,O2,t) ^ t [(3t(Oj)-xtaJt)) '|z^l 

and prove two crucial lemruis. 
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Lemna 6{ The joint conditional density 


,x (Oj) ^ ^ ^ is G&ussian with lionrandom 

conditional cross-covariance 


Prwf : First, the conditional density is 

Gaussian because and are jointly Gaussian 
random processes. Assume ^ ^ 2 * then 


tdiere » (zfs), 0^ < s < t>. 


Each of the densities in (13) is the result of a 
linear smoothing operation; hence, each is 
Gaussian with lionrandom covariance 

Jo PCG' 2 ,a 2 ,t) , respectively 120). 

X 2 


Thus the cross^covariance aatiafles 

**0 lo 
X 2 

•p"^(02,a2,t) 

and PCO.jOj/t) is also nonrandom (here P(0 ,0 ,t) 
is invertible because is stabllizable) , 

Lemma 6 allows the off-line ccanputation of . 

of P(g^»<^ 2 '*'^ equations of KWakemiaak [24] 

(forOj^^o^); 

P(o^,02-t) » P(o^) ¥'{a^, 0^) 
t 

-P(O^) [ j ¥'(T, 0^) H' (T) r"^(T). 

. H(T) I* (T, 02>dT]P(a2) 
(14) 


'l'(t) - [F(t.) - P(t)H'(t) R"^(t)H(t)] ’I' (t); 

v (0) - I as) 

%i4iere the Kalman filter error covariance matrix 
P(t) “ P(t, t, t) is computed by the Riccati 
equation 


P(tJ " P(t) P(t) + P(t) P’(t) + Q(t) 

-p(t)H'(t) R*^(t) H(t)p.(t) as) 

r . emiaa 7 ; The conditional cross-covariance 
satisfies 

P(0, t, t) « K(t, O) P(t) (17) 

where 

|^K'(t,a) - -lF'(t) + p'^(t) Q(t)] K'(t, a); 

KMo, o) = I (18) 


Proof ; Let 

P(0,t) " E[(x(a) - 5(a|ci)) Cx(t) - x(t|t))*], and 
consider 

P(o,t,t) - p(0,t) E'^((iJ(ojcri - x(o|t)) 

(x(t) - X(tjt)) •] (19) 


Since x(o]o) - S(olt> is measurable with respect 

t 

to z , the projection theorem implies that (19) 
equals zero. The proof is conclude by noticing 
that P(c,t) « K(t,a) p(t) [21]. 


Returning to the proof of Theorem 1 for k«3, we 
first augment the state of (2) with 
y^^(t), i«l,...,4 (if a particuleu: y^ can be ob- 
tained as a linear combination of the xls, we 
need not augi ent the state with that y^) . Then 
the Kalman F: Iter for the system (2), (4), ( 11 ) 
generates x(t|t) and y(t|t). We define the 8x8 
conditional covariance matrix 


V(0j^,02.t) 


|'p(V2,t) 





wtu»re 


s(aj,02,.) 
T(0^,0^,l ) 





SCOj.O^.t) = E^[(x(0^) - S(0^|t)). 

.(y(02) - y.CO^lt))’! 


P(aj^,T,t)dT 


T(0j.02,t) = E^f(y(a^} y(Oj_!t)) . 

*^2 . 


( 20 ) 


= 1 ] 
0 0 


P(T, ,T2,t)d-r^dT2 


( 21 ) 


v;e also dc-fine T(t) « T(t.t,t), S(t) = s(t,t,t), 
V(t) « V(t,t,t). 

The (characteristic function of a Gaussian 
random vec tor x vjith ncan n and covariance P is 
given by 


Jd^(u) « Ji;[exp (iu'x) ] «* expfi u’n - j u'pul (22) 


Thus (for j“l,2,3) 


It can alj-o be alio\/n that 


y^(tlt) + 


(23) 


X^2(t|t) » (YjCtlt) + 


y^(t|t) + 2 Tii(t) 


(24) 

y^(tlt) .+ I T^^(t) 


(25) 


X23(t|.t) = (y^(t|t) + 


E^(e y 3 (t)) = ly^ctjt) + T^ 3 (t)) . 

. e 

(2G) 

since (23)- (20) represent instantaneous nonlinear 
functions of 7 (tlt)> they can be computed v;ith a 
finite dii jnsional ostinator. 


of the form 

W(t,0) = U(t)V(tJ) 

then a(t) is the output of a finite ciiinonsional 
linear system driven by x(t) tl9], and a(t|t) 
can be computed by the linear finite dimensional 
Kalmern filter for the augmented state 
(x'(t) a'(t))'. 

It can easily be shown that • 
t 

E [e ■' ^^(o^) x^(02)) 


(28) 


t|t) + ^ T^^ct) 


[p^4(o^,a2,t) 


(29) 

Hence 

E^'te ^ Kit)] 

. + T,- r^^(t) 


. (E 


+ B 


: r t ,«^1 

_t 


X2(0^) X^(02) da^dCJ^ 


0 0 

t Oi' ■■ 


Si2<t,o^,t) x^(02) da^ 


LOO 


+ E 


a ‘ 


If 

Loo 


t o 




I Sj^otfc.Oj^'t) (30) 

0 0 

Consider the :.econd tem in (30) . Lcnr.a 7 implies 
that V(o,t,t) can he v/ritten as 


Uov; consider 5(t) (see (12)) j the approach 
here will 'ir: the rorhiction of this problon to 


a(tlt) = E"[a(t>). = 


\7(t,0) ::{O)d0 


(27) 


V(0,t,t) = L(t,0) V(t) 


where 




• •(u 4 j, 

• ’I. ' (M ''(h) 


1/ a.O) ; 
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..-1 , A 

'O^l(t) 


V -^(t) = U(t) * 

-U^2'(t) 

U22(t). 

Since Sj ^2 is the 

(1/6) element of 

that 




t 


s^2(t,ai.t) da^da^ 


Lo 0 


I y^(o,) da^j 

t 1 

I C^ L(t,o) y^j(a)daJ v(t)e^ 


= E^{6' (t)] v(t)e. 


th 


(32) 


(where is the j unit vector) can be 
computed v^ith a finite dinensionai linear 
estinvator. 

Similar reasoning implies that the third term 
in (30) 




Ji 

1*0 0 




r ty ^1 




0 


®s ^<t,a2)d02 I X^(a^)d0^ 


• V(t)e^ = E^te' (t)] V(t)e^ 


(33) 


Can be corputed v/ith a finite dimensional linear 
estimator. 


Let Y be defined by 

dY(t) = x.,(t) y^(t) dt; Y(0) « 0 


(3d) 


Then the nonlinear filtering equation 122] yields 
dY(t[t) = L'^(x,(t> y^(t)] dt 

+ {E*^[Y(t) x'(t)J - Y(t|t) x'(t[t)}. 

■H* (t) r. ^(t) [dz(t)-II(t) x(tlt)dt] (35) 

■Die first term in (30) (the HriTt tenn) can be 
written as 


E^tx 2 (t) y^(t)] => X 2 (t|t) y^(t|t) + S 2 ^(t) 

v/hich is computable "finite dimensionally". 
We have that the gain term 

E^lY(t) x’(t)J - Y(tjt) x'(tjt) 

t 


( 36 ) 


0 0 


- x' (tlt)E [X2(0j^)x^(02)] j da^dOj^ 

■r 




0 0 


P(t) 


j K y^(0^)da^ 

^0 J 

t 

^ [1 { K^(t.°2^ <302^^1 


(E^teMt.) + E^CO' (t)]}P(t)' 


(37) 


where K. denotes the i— row of K. since 

P(0^«tft) and J P( 02 »t/t) are both separable! 

0 

(37) and hence (35) can be computed finite di^n- 
slonally. The proof for Jc«3 is now complete* 

The optimal isstimator for the 3x3 nili>otent system 
(10) consist'! of a Kalman filter for the augmented 
state consisiing of x, y, o, 8, 6, and e 
(defined in ;2), (11), (37), (37), (32), and (33) 
respectively with observations (4) , together with 
the nonllnea; stochastic equation (35) . These 
are followed by nonlinear transforiTiations as shown 
in (23)-(26) and (30). 

A block diagram illustrating the nonlinear esti- 
mator is shov.n in Figure 1; notice that the esti- 
mator Is driven by the innovations process 
dv(t) - dz(t) - H(t) xCtjt) dt. 


5. rji:r:M'Ai,j/.Arioii!» 


'r»K» r*. I !•»%*!». 

I < O -1 (J f I , 


OF tse 

IS POOR 



Theorem 2 ; Consider (l}-(4), where ^ is 


solvable and is nilpotent. Then the condi- 
tional mean X(tjt) can be computed with a finite 
dinensionaJ system of nonlinear stochastic 
differential equations. 

If, however, ^is solvable and is not 
nilpotent, the optimal estimator will be Infinite 
dimensional. For example, if k«2 we must compute 
the conditional expectation 



(see [16) for further details) . 


A proof similar to that for Theorem 2 yields a 
further qenaraliaation. 

Theorem 3 : Consider the equations (2)- (4) 

and 

N 

X(t) (A (t) +• J u (t)A.)X(t)j X( 0 ) - 1 ( 38 ) 

” i -1 ^ ^ 

Let^ - Ajj(t) (Vt)}^, and let 

be the ideal in generated by {a^, . . . ,Ajj}, 

Assume that is nilpotent. llien the condi- 

tional mean X(t|t) can computed with a finite 
dimensional system of nonlinear stochastic 
differential equations. 

Remark : Notice th it if A^Ct) is time-varying, 

the nilpotency of does not imply that X(t) 
evolves on a solvable Lie group. 

Theorem 1 can also be extended to other systems 
with polynomial nonlinearities (the proof is 
similar) 

Theorem 4 ; Consider the linear system des- 
cribed by (2) and (4), and define 


t 0^ 1 

Y(t) “ 1 I ••• 1 ) ••• >• 

J J i ^1 "l “i 

0 0 0 


• "A > ••• 


Where are arbitrary deterministic matrix - 

valued functions. The subscripts 

{n } are not necessarily distinct, and i and I 
are not necessarily less than or equal to k. 

r 1 

Then and E^^Le Y(t)J can be conqputed 

with a finite dimensional system of nonlinear 
stochastic differential equations. 

6, RATH- INDEPENDENT INTEGRALS 


The results of this section are based upon the 
work of Brockett [19] au^d Gruber [13] on stability 
analysis using exact differentials and path-inde- 
pendent Integrals. We will consider equation (2) 
with x(t) Cr", in . which it will be assumed that 
each component x^(t) is - tiroes mean-square 
(m.s) differentiable (m^ nay be equal to zero for 
some 1)} let k., be the number of components of x 
%daich are at l<^ast j-titnes m.s. differentiable.' 
Also, we deflni.! to be the kj x n selector 
matrix vdilch selects the coa^nents of x that are 
at least j-times m.s. differentiable! 

®i,l " 1 1 J ^ 3 


' » 1 if m. > j and m. < j for i ,+l<a<i -1, 

>,i i - f . p-1 p 


for p>l 


S =0 otherwise 


Finally, let 


5j(j) d 


4 [s^x,.,] 

dt-* 


( 39 ) 


If n “ max (m. ) and f is a continuous function, 
Ki<n ^ 


the random process 







is said to be independent of path (in the mean- 
square) if there exists a function g such that 


Y(t) 


g(x(t) 


.(n-1) 


(t),x(0) 




( 0 )) 

(41) 


where equality in (41) is in the mean— square 
sense (these definitions could also be placed in 
the "almost sure" framevrark) . 


First we consider the case n=l. If g is twice 
continuously differentiable, and x(0) is kno\vn, 
then Y(t|t) satisfies the nonlinear filtering 
equation [ 22 J 
'' . 

dY(t|t) « {e [g^(x(t))F(t)x(t)] 

+ i E^(tr(Q(t)g (x(t)))]}dt 

« ywC 

+ {E^(Yft)x' ft) I - Y(t|t)x' (tft)} . 

•H* (t)R ^(t) {dz(t)-II(t)x(t|t)dt] (42) 
where g is the gradient and g is the matrix 

X XX 

of second partials. Since x and its deriva- 
tives are Gaussian, it is easy to see that (42) 
be computed in .terms of the conditional mean 
and covariance of x; thus the estinvate Y{t|t) of 
a path-independent Integral can be computed with 
a finite dimensional nonlinear estimator. This 
result can obviously be extended to the ease in 
which m > 1. 

Bxample 3 ; Define y^ as in (11). If 
' t 

Y(t) » Ix. (o)y. (a) + X.. (o)y. (o)lda 
} ^ J 31 - 

0 

^ (t) -y^(0)y^(0) 

then Y(tjt) is finite dimensionally computable.. 

Exangle 4 i Assume x (0) « ^0, and let 
t 

Y(t) = j < T(a)s^x(o) + T(a>^^^\(o), 

° 1 

T(a)s x{a) > do 


Theorem >i Consider the linear system descri- 
bed by (2) .ind (4), where x(0) is known. Let 
B(t) • Y(t) + 6(t) (43) 

where Y is i patli- independent integral defined by 
(40), and E''(6(t)] = 0. Then Y(t|t) can be com- 
puted with a finite dimensional system of nonlinear 
stochastic differential equations. 

Finally, we state a theorem concerning the inte- 
grals of qut-dratic differentials, which is based 
on a result of Gruber' [13}. 

Theorem C ; Consider the linear system described 


by (2) and ( 4 ), where x(0) is known. Define the 
scalar differential operator 


P(D) 



n-1 




(44) 


, , m r i 

p(s) » s + i PqS (45) 

i“0 

and the x k^ matrix, differential operator 

Q(D) « ? Gi ^ (46) 

i"0 dt 


Q(s> « I Q.s^ (47) 

i-O ^ 

where q < m. Assume that the matrix R(s) satisfies 

P(s) Q' (s) + p(-s) 0 (s) s R' (-S) S(s) ( 48 ) 

Then 

Y(t) =1 2 < p(D)s”'x(0), Q(D)s"x(0) > 

° - < R(D)s”x(0), R(D)s’^x(0) > da 


is independont of path, and ■y(t[t) can be computed 
with a finite dimensional system of nonlinear 
stochastic differential equations. 


where denotes inner product. Then 

Y(.t) < T(t)S^x(t), T(t)S^x(t) > 

and Y(tlt) is finite dincnoionally computable. 

A simple extension of these ideas is the foll'jwing. 
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